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We compare the dark energy model constraints obtained by using recent standard ruler data 
(Baryon Acoustic Oscillations (BAO) at z = 0.2 and z = 0.35 and Cosmic Microwave Background 
(CMB) shift parameters R and l a ) with the corresponding constraints obtained by using recent Type 
la Supernovae (Snla) standard candle data (ESSENCE+SNLS+HST from |astro-ph/0701510| . We 
find that, even though both classes of data are consistent with ACDM at the 2a level, there is a 
systematic difference between the two classes of data. In particular, we find that for practically 
all values of the parameters (f2om,fit) in the 2a range of the the 3-year WMAP data (WMAP3) 
best fit, ACDM is significantly more consistent with the Snla data than with the CMB+BAO data. 
For example for (S7om, fi(>) = (0.24, 0.042) corresponding to the best fit values of WMAP3, the dark 
energy equation of state parametrization w(z) = wo + wi best fit is at a 0.5a distance from 
ACDM (w = -l.ioi = 0) using the Snla data and 1.7a away from ACDM using the CMB+BAO 
data. There is a similar trend in the earlier data (SNLS vs CMB+BAO at z = 0.35). This 
trend is such that the standard ruler CMB+BAO data show a mild preference for crossing of the 
phantom divide line to = — 1, while the recent Snla data favor ACDM. Despite of this mild difference 
in trends, we find no statistically significant evidence for violation of the cosmic distance duality 
relation n = dA ( z ^i+ z yt — 1- For example, using a prior of Qom = 0.24, we find n — 0.95 ± 0.025 in 
the redshift range < z < 2, which is consistent with distance duality at the 2a level. 

PACS numbers: 



I. INTRODUCTION 

The accelerated expansion of the universe has been 
confirmed during the last decade by several observational 
probes 0, 0, S, H, i, H, . The origin of this acceleration 
may be attributed to either dark energy with negative 
pressure, or to a modification of General relativity that 
makes gravity repulsive at recent times on cosmological 
scales. In order to distinguish between these two possi- 
bilities and identify in detail the gravitational properties 
of dark energy or modified gravity two developments are 
required @: 

1. Detailed observation of linear cosmic density per- 
turbations 5(z) = Sp ^ z ' at recent redshifts. 

2. Detailed mapping of the expansion rate H(z) as a 
function of the redshift z. 

The later is equivalent to identifying the function w(z) 
defined as 

I dMSHjzf) 
w(z) = -l+-(l + z) dlnz , (1.1) 

where 8H(z) 2 = H{z) 2 /H% - £7 0m (l + zf - n 0r (l + z) 4 
accounts for all terms in the Friedmann equation not re- 
lated to matter and radiation. If the origin of the ac- 
celerating expansion is dark energy then w(z) may be 



identified with the dark energy equation of state param- 
eter w(z) = The cosmological constant (w(z) = —1) 
corresponds to a constant dark energy density. 

It has been shown [9( that a w(z) observed to cross 
the line w(z) = —1 (phantom divide line) is very hard 
to accommodate in a consistent theory in the context of 
General Relativity. On the other hand, such a crossing 
can be easily accommodated in the context of extensions 
of General Relativity 0]. Therefore, the crossing of the 
phantom divide line w = — 1 could be interpreted as a 
hint in the direction of modified gravity. Such a hint 
would clearly need to be verified by observations of linear 
density perturbation evolution through eg weak lensing 
[To| or the redshift distortion factor [ill ]. 

There are two classes of probes that may be used to 
observe the expansion rate H(z) or equivalently w(z) 

• Standard Candles are luminous sources of known 
intrinsic luminosity which may be used to measure 
the luminosity distance which, assuming flatness, is 
connected to H(z) as 

d L (z)=c(l + z) J Z dz'j^- (1.2) 

Useful standard candles in cosmology are Type la 
supernovae P, [H, HI, [I3l (Snla) and the less accu- 
rate but more luminous Gamma Ray Bursts [l5j 
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• Standard Rulers are objects of known comoving size 
which may be used to measure the angular diameter 



2 



distance which, in a flat universe, is related to H(z) 
as 



d A (z) 



I ■ z 



(1.3) 



The most useful standard ruler in cosmology is 
the last scattering horizon, the scale of which can 
be measured either directly at z ~ 1089 through 
the CMB temperature power spectrum or indi- 
rectly through Baryon Acoustic Oscillations (BAO) 
on the matter power spectrum at low redshifts. 
Clusters of galaxies [H, QjJ and radio galaxies fl8j 
may also be used as standard rulers under cer- 
tain assumptions but they are less accurate than 
CMB+BAO. 

Early Snla data put together with more recent such 
data through the Gold dataset [l], [l2j have been used 
to reconstruct w(z), and have demonstrated a mild pref- 
erence for a w(z) that crossed the phantom divide line 
(l9l I20I ]. A cosmological constant remained consistent 
but only at the 2a level. However, the Gold dataset 
has been shown to suffer from systematics due to the in- 
homogeneous origin of the data [2l|. More recent Snla 
data (SNLS[I3, ESSENCE^, HST[T3]) re-compiled in 
[l3l ] have demonstrated a higher level of consistency with 
ACDM and showed no trend for a redshift dependent 
equation of state. 

On the other hand, the use of standard rulers 
(CMB+BAO) has rarely been studied independent of 
Snla due to the small number of datapoints involved (see 
however [H, [H [H|). It has been pointed out [II] that 
the latest BAO data "require slightly stronger cosmo- 
logical acceleration at low redshifts than ACDM " . This 
statement is equivalent to a trend towards a w{z) < — 1 
at low z, and therefore a possibility of crossing the PDL 
w = —1. The goal of this paper is to quantify this state- 
ment in some detail by comparing the best fit form of 
w{z) obtained from the Snla data to the corresponding 
form obtained from the CMB+BAO data. This compari- 
son is done quantitatively by identifying the quality of fit 
of ACDM in the context of each dataset. In particular, 
we consider the Chevalier- Polarski-Linder (CPL) [27L [28l] 
paramctrization 



w(z) = wq + w\ 



l + z 



(1.4) 



and, assuming flatness, we identify the "distance" in 
units of a (cr-distance) of the parameter space point 
(wo,wi) = (—1,0) corresponding to ACDM from the best 
fit point (wq, wi) for each dataset (Snla standard candles 
or CMB+BAO standard rulers) and for several priors of 
(Oom,^b)- We thus identify an interesting systematic 
difference in trends between the two datasets. 

We also discuss the implications of this difference in 
trends on the distance duality relation 



r)(z) 



d L (z) 



d A {z)(l + zf 



1 



(1.5) 



which measures quantitatively the agreement between lu- 
minosity and angular diameter distances. This relation 
has been shown to be respected when clusters of galaxies 
are used as standard rulers 1171. 



II. LIKELIHOOD CALCULATIONS 

We assume a CPL parametrization for w(z) and apply 
the maximum likelihood method separately for standard 
rulers (CMB+BAO) and standard candles (Snla) assum- 
ing flatness. The corresponding late time form of H(z) 
for the CPL parametrization is 

H 2 (z) = H<-[n 0m (l + zf + 

+ (l-f] nm )(l + ^) 3(1+ ™ 0+ffil) ^] (2.1) 

At earlier times this needs to be generalized taking into 
account radiation ie 

E 2 {a) = = n m ( a + a eq )a~ 4 + fl de X{a) (2.2) 

where a = 1/(1 + z), fide = 1 — &m — ^rad and 



X(a) 



.3 r£±-(«2W 



exp 



a ~ 3(l+wo-t-«>i) e — 3«;i(l— a) 



(2.3) 



with the CPL parametrization w(a) = wq + Wi(X — a). 



Standard Rulers 



1. CMB 



We use the datapoints (R, l a , Sl^h 2 ) of Ref. [25[ where 
R, I a are two shift parameters: 



The scaled distance to recombination 



m 



R = \ ^0m— Y r(zcMB) 



(2.4) 



where t(zcmb) is the comoving distance from the 
observer to redshift z and is given by 



r(z) 



(2.5) 



dz 

HvJo W) 
with E(z) = H(z)/H . 

• The angular scale of the sound horizon at recombi- 
nation 



_ r^acMB) 
r s (acMB) 



(2.6) 



where r s (acMB) is the comoving sound horizon at 
recombination given by 



c f 
-"0 Jo 



c s {a) 
a 2 E(a) 



da (2.7) 
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with the sound speed being c s (a) = 1/^/3(1 + R a) 

and acMB — j-rr 1 , where z C mb = 1089. Actu- 

ally , zcmb has a weak dependence on il m and f2f, 
[23 ] but we have checked that the sound horizon 
changes only to less than 0.1%. The quantity R^, 
is actually the photon-baryon energy-density ratio, 
and its value can be calculated using R = | = 
31500fi 6 /i 2 (T CA/B /2.7is:)- 4 . 

For a flat prior, the 3-year WMAP data (WMAP3) Q 
measured best fit values are [25| 

R \ / 1.70 ±0.03 
' c 'MB = | L = 302.2 ±1.2 I (2. 
n b h 2 J \ 0.022 ± 0.00082 

The corresponding normalized covariance matrix is [2£ 



pnorm 
^CMB 



1 -0.09047 -0.01970 
-0.09047 1 -0.6283 | (2.9) 
-0.01970 -0.6283 1 



from which the covariance matrix can be found to be: 

(CcMB)ij = {CcMB)ij <J v^ MB a v^ MB ( 2 - 10 ) 
where a-wi are the la errors of the measured best fit 

v CMB 



values of eq. 

We thus use equations (f2T8|) . d2~4]) and (|2~6]) to define 



where the dilation scale 
Dv(zbao) = 



dz 



zbao 



H(z)J H(z BAO ) 



1/3 



(2.15) 



encodes the visual distortion of a spherical object due to 
the non-Euclidianity of a FRW spacetime, and is equiv- 
alent to the geometric mean of the distortion along the 
line of sight and two orthogonal directions. We thus con- 
struct 



Xbao = 



r a (Zdcc) 

D v (0.2) 

r s(Zdcc) 

D v (0.35) 



0.1980 
- 0.1094 



and using the inverse covariance matrix [26 



'BAO 



35059 -24031 
-24031 108300 



we find the contribution of BAO to x 2 as 



Xbao 



Xbao T Cbao 1 Xbao 



B. Standard Candles 



1. Snla 



(2.16) 



(2.17) 



(2.18) 



R- 1.70 

Xcmb = j l a — 302.2 
n b h 2 - 0.022 



(2.11) 



and construct the contribution of CMB to the % 2 as 

Xcmb = ^cmb T Ccmb ^cmb (2-12) 

with 

/ 1131.32 4.8061 5234.42 
Comb -1 = 4.8061 1.1678 1077.22 )(2.13) 
\ 5234.42 1077.22 2.48145 x 10 6 

Notice that Xcmb depends on four parameters (fionn 
Wq and Wi). Due to the large number of parameters 
involved, in what follows we will consider various different 
priors on the parameters f^om, 

2. BAO 

As in the case of the CMB, we appl y t he maximum 
likelihood method using the datapoints [26| 



( r D ZC (0 2) } = 0-1980 ± 0.0058 \ 
VBA ° = (#^ = 0.1094 ±0.0033 J ' (2 ' 14) 



We use the Snla dataset of Davis et. al. 13] consisting 
of four subsets: ESSENCE [H (60 pomtsVSNLS Q (57 
points), nearby [l[ (45 points) and HST ( 30 points). 

These observations provide the apparent magnitude 
m(z) of the supernovae at peak brightness after imple- 
menting the correction for galactic extinction, the K- 
correction and the light curve width-luminosity correc- 
tion. The resulting apparent magnitude m{z) is related 
to the luminosity distance Dl(z) through 

m th {z) = M(M, H ) + 5log w (D L (z)) (2.19) 

where in a flat cosmological model 



D L {z) = (1 + z) / dz' 



Hn 



o i/(z';ri m,w ,wi) 



(2.20) 



is the Hubble free luminosity distance (Hod^), and M is 
the magnitude zero point offset and depends on the abso- 
lute magnitude M and on the present Hubble parameter 
H as 



M = M + 5Z <?io(|M + 25 = 
Mpc 

= M — blog w h + 42.38. 



(2.21) 



The parameter M is the absolute magnitude which is 
assumed to be constant after the above mentioned cor- 
rections have been implemented in m(z). 
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FIG. 1: The 68.3% and 95.4% \ 2 confidence contours in the (wo,^i) parameter space for each dataset category: Fig la for 
Snla (based on data from 43]) and Fig lb for CMB+BAO data (based on data from [H[H) for n <>m = 0.24 and fi 6 = 0.042 
(best fits of WMAP3). The blue dots correspond to the (wo,u>i) best fit, while the yellow dots correspond to ACDM (-1,0). 



The Snla datapoints are given, after the corrections 
have been implemented, in terms of the distance modulus 



Hobs(zi) = m obs (zi) - M 



(2.22) 



The theoretical model parameters are determined by 
minimizing the quantity 



2 / \ (fJ-obs(zi) - ^thiZi)) 2 . . 

XSn/a(^0m,W0,Wl) = 2^ 2 ( 2 ' 23 ) 

i=l °M i 

where N — 192 and a 2 , „■ are the errors due to flux uncer- 
tainties, intrinsic dispersion of Snla absolute magnitude 
and peculiar velocity dispersion. These errors are as- 
sumed to be Gaussian and uncorrelated. The theoretical 
distance modulus is defined as 

Hth(zi) = m th (zi) - M = 5log w {D L (z)) + ^ (2.24) 

where 

Mo = 42.38 - Uog w h (2.25) 



and [lohs is given by (|2.22|) . The steps we followed for the 
minimization of (|2.23p are described in detail in Refs. 



C. Results 

We consider separately the standard ruler data (xsr = 
Xcmb + Xbao) and the standard candle data (xsnia)> 
and perform minimization of the corresponding \ 2 with 



respect to the parameters Wq and w\ for various priors 
of f^om and in the 2a range of the the WMAP3 best 
fit ie 0.17 < fiom < 0.31, 0.034 < fl b < 0.049. 

In Fig. 1 we show the 68.3% and 95.4% x 2 confidence 
contours in the (wq , w\ ) parameter space for the two 
dataset categories (standard ruler and standard candle 
data) for Q 0m = 0.24 and Q b = 0.042 (the best fit of 
the WMAP3 CMB data 0). Fig. la shows the (w ,w x ) 
contours obtained using Snla data [13| (standard candles) 
while Fig. lb shows the corresponding contours assum- 
ing CMB+BAO data [Hill (standard rulers). The blue 
dots correspond to the (wo,wi) best fit, while the yel- 
low dots correspond to ACDM (wq,wi) = (—1,0). The 
distance in units of a (er-distance d a ) of the best fit to 
ACDM was found by converting Ax 2 = Xacdm ~ xLm 
to d a ie solving [3l| 



l-r(l,A X 2 /2)/r(l) = Erf(d c 



(2.26) 



for d a (cr-distance) , where A^ 2 is the \ 2 difference be- 
tween the best-fit and ACDM and Erf () is the error func- 
tion. Notice that ACDM is consistent at less than la 
level according to the Snla data (d^ nIa ~ 0.5 in Fig. 
la), while the corresponding consistency level reduces 
to df R ~ 1.7a for the standard ruler CMB+BAO data. 
This mild difference in trends between standard candles 
and standard rulers persists also for all values of fiom in 
the 2a range of WMAP3 best fit. This is demonstrated 
in Fig. 2 where we show the cr-distance df R superposed 
with d s a nIa as a function of n 0m for Q b = 0.034 (Fig. 
2a), fl b = 0.042 (Fig. 2b) and Q b = 0.049 (Fig. 2c). 
These values of tt b span the 2a range of the correspond- 
ing WMAP3 best fit. Notice that the cr-distance between 
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FIG. 2: The a-distance of the best fit parameters (wo,Wi) to ACDM (-1,0) for Snla data (black lines) and for CMB-BAO data 
(colored lines), as a function of fi 0m for = 0.034 Fig. 2a, Q b = 0.042 Fig. 2b and Q b = 0.049 Fig. 2c. 
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a= 0.0344 
A = 0.0382 
fl b = 0.042 
fi b = 0.0458 
fi b = 0.0496 



FIG. 3: The ratios of the <r-distances of the best fit parameters (ti?o,wi) to ACDM (-1,0) for Snla standard candles over the 
one for CMB-BAO standard ruler data as defined in eq. (|2.27[) , as a function of ilom for various values of fib- 



best fit values and ACDM values is consistently larger 
when using standard ruler data (colored lines are consis- 
tently above black lines). 

An alternative way to see this trend is to plot the ratios 
r of the cr-distances defined as: 

Kflom) - J^r (2-27) 

These plots are shown in Fig. 3 for five values of Of, span- 
ning the 2a range of WMAP3. Notice that the colored 
lines are consistently above the line r = 1 indicating that 
the er-distance is found to be consistently larger when us- 
ing standard ruler data. An exception to this rule is the 



case corresponding to high values of both Sit, and f2om set 
to values 2a or more, away from their best fit (see ma- 
genta line corresponding to Qb = 0.049, for Q 0m > 0.29). 

The above plots reveal a consistent trend of the stan- 
dard ruler CMB+BAO data for a mild preference for 
crossing of the phantom divide line w = — 1, while the re- 
cent Snla data seem to favor ACDM. Figs. 4 and 5 show 
corresponding plots obtained using earlier data (SNLS 
vs CMB+BAO H,i|atz = 0.35), where a similar 
consistent trend is observed. 

An interesting feature of the contours of Fig. lb is the 
deformation appearing for relatively large values of w\ . 
There is a simple way to understand this deformation. 
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FIG. 4: The 68.3% and 95.4% x 2 confidence contours in the (u>o,wi) parameter space for each of the old dataset categories [l4l. l33l] 
for fiom = 0.24. The green dots correspond to the (wo,Wi) best fit, while the yellow dots correspond to ACDM . 




FIG. 5: (a) The ratio of the er-distance of the best fit parameters (wo, Wi ) to ACDM (-1,0) for Snla over the one for CMB-BAO 
data, as a function of fio m for the early CMB+BAO data[33l]. (b) The a-distance of the best fit parameters (u>o,Wi) to ACDM 
for Snla data (black line) and for CMB+BAO data (red line), as a function of fiom- 



For (u>o,wi) parameter values satisfying u>o + w± ~ 0, 
the dark energy equation of state is approximately con- 
stant and positive at early times and therefore the corre- 
sponding dark energy density dominates over the matter 
density. As a result the Hubble expansion rate is sig- 
nificantly modified over the whole range from z = to 
zcmb and the corresponding integral of the shift param- 
eters becomes very sensitive to parameter changes. The 
effect is even more significant for the shift parameter l a 
which involves the sound horizon r s in the denominator 
(see eq. (I2.6[) . The sound horizon drops more dramati- 
cally than the shift parameter R when the dark energy 
dominates at early times because the corresponding in- 
tegral (|2.7p depends only on the early time behavior of 
the expansion rate H(a). This effect is demonstrated by 
plotting the wo + w\ = line in Fig. lb which coin- 
cides approximately with the region where the contour 



deformation starts. The same line is also plotted in Fig. 
4b corresponding to early CMB+BAO dataJHJ and in- 
volving only one shift parameter (R). In this case the 
deformation effect is milder because the z integral cor- 
responding to R spreads over a wide range of redshifts 
from z = to z = 1089 and the effect of dark energy 
domination is somewhat smeared out. 



We also construct the likelihood contours using the 
combined Snla+CMB+BAO data for f2 0m = 0.24 and 
fib = 0.042 corresponding to the best fit WMAP3 param- 
eter values. As expected, the cr-distance between best fit 
and ACDM is at about la, ie intermediate between the 
standard candle and standard ruler cases (see Fig. 6). 
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FIG. 6: The 68.3% and 95.4% % 2 confidence contours in 
the (wo,wi) parameter space for the combined datasets 
Snla+CMB+BAO for Sl 0m = 0.24. The blue dot corresponds 
to the (mi ,wi) best fit while the red dot to ACDM (-1,0). 

III. CONCLUSIONS-DISCUSSION 

We have demonstrated that there is a systematic dif- 
ference in trends between standard candle (Snla) and 
standard ruler (CMB+BAO) data. The former data are 
significantly more consistent with ACDM than the later 
for practically all (f2om,^fi) parameter priors within the 
2<T range of WMAP3. In fact, the standard ruler data 
demonstrate a mild preference for a best fit w(z) that 
crosses the phantom divide line w = — 1. 

This systematic difference in trends can be attributed 
to one of the following: 

• Statistical Effects: There is an (l]j] m , param- 
eter range where both datasets are consistent with 
each other and with ACDM at the 2a level (see eg. 
Fig. 3b with Slo m ~ 0.25). Therefore, for these pa- 
rameter values the two datasets are consistent with 
each other and with ACDM at the la — 2a level 
and the trend we observe could well be a statistical 
fluctuation. 

• Systematic-Physical Effects: As discussed in 
Rcf. 34j] distances based on standard candles and 
standard rulers should agree as long as three con- 
ditions are met: (1) photon number is conserved, 
(2) gravity is described by a metric theory and (3) 
photons are traveling on unique null geodesies. If 
at least one of these conditions is not met then 
equations (|1.2[) and (|1.3|) will lead to generically 
different forms for the Hubble expansion rate H{z) 
due to the violation of the distance duality relation 



(jl.5[) . For example, lensing of Snla by compact ob- 
jects, if not properly accounted for, would tend to 
violate condition (3) and induce artificial brighten- 
ing of distant Snla. Alternatively, photon number 
violation (due eg to photon mixing [351 ]) would lead 
to artificial dimming of the Snla. 
In order to investigate the possible existence of system- 
atic physical effects, we have used our results to test the 
cosmic distance duality relation (|1.5|) . In particular, we 
use our results for the best fit parameter values ((wq, w±)) 
and their error bars obtained from each dataset to derive 
constraints on the parameter rj(z). These constraints 
are shown in Fig. 7a for flo m = 0.24 and in Fig. 7b 
flom = 0.27. Clearly, the anticipated value rj = 1 is 
within 2a for both priors used. Assuming a prior of 
flom = 0.24 and taking an average value for t)(z) in the 
range < z < 40 (as for large enough z r\(z) converges 
(see Fig 7)), yielding the value f) — 0.96 ± 0.07 which is 
within la from the anticipated value 77 = 1. The con- 
sistency is somewhat reduced if we average over a more 
recent redshift range. In the range 1 < z < 2 we find 
fj = 0.95 ± 0.025 which is consistent with the anticipated 
value rj = 1 at the 2a level. Similar results are obtained 
for other priors of f2o m within 2a from the WMAP3 best 
fit. Therefore, despite the mild difference in trends be- 
tween Snla standard candles and CMB+BAO standard 
rulers, we find no statistically significant evidence for vi- 
olation of the distance duality relation. 

An interesting extension of this work would be the in- 
clusion of more data from both categories. For example 
gamma ray bursts [l5| could also be included as stan- 
dard candles and X-ray profiles of clusters [l7| or radio 
galaxies [l8| could be included as standard rulers in or- 
der to investigate if the mild difference in trends we have 
identified, persists in more general categories of data. 

The Mathematica files with the numeri- 
cal analysis of the paper can be found at 
http://leandros.physics.uoi.gr/rulcand/rulcand.htm 



Acknowledgements 

We thank Y. Wang, P. Mukherjee, W. Percival and 
E. Majerotto for useful discussions. This work was sup- 
ported by the European Research and Training Network 
MRTPN-CT-2006 035863-1 (UniverseNet), by the Uni- 
versity of the Basque Country through research grant 
GIU06/37, and by the Spanish Ministry of Education 
and Culture through research grant FIS2004-01626. S.N. 
acknowledges support from the Greek State Scholarships 
Foundation (I.K.Y.). 



[1] A. G. Riess et al. [Supernova Search Team Collabora- [2] D. N. Spergel et al. [WMAP Collaboration], Astrophys. 
tion], Astrophys. J. 607, 665 (2004). J. Suppl. 170, 377 (2007) [arXiv:astro-ph/0603449] . 



8 




[3] 
[4] 
[5] 
[6] 

[7] 
[8] 



[9] 
[10] 



A. C. S. Readhe ad et al, Astrophys. J. 609, 498 (2004) 
[arXiv: |astro-ph/0402359l . 

J. H. Goldstein et al, Astrophys. J. 599, 773 (2003) 
arXiv:astroph/0212517 . 

Rebolo et al., MNRAS 353 747R (2004) 



arXiv:astro-ph/0402466 



M. Tegmark et al. [SDSS Collaboration], Phys. Rev. D 
69, 103501 (2004) [arXiv:astro-ph/0310723l 
E. Ha wkins et al, Mon. Not. R oy. Astron. Soc. 346, 78 
(2003) [arXiv:astro-ph/0212375] 

B. Boisseau, G. Esposito-Farese, D. Polarski and 
A. A. Starobinsky, P hys. Rev. Lett. 85, 2236 (2000) 
arXiv:gr-qc/0001066 ; E. Bertschinger, Astrophys. J. 
648, 797 (2006) |arXiv:astro- ph/0604485 ]; J. P. Uzan , 
Gen. Rel. Grav. 39, 307 (2007) farXiv:astro-ph /0605313| ; 
J. P. Uzan and F. Bernardeau, Phys. Rev. D 64, 
083004 (2001) [arXiv:hep-ph/00 12011 ; L. Perivolaropou - 



arXiv:astro-ph/0504582 ; 
A~. Melchiorri, Phys. 



los, JCAP 0510, 001 (2005) 
R. Caldwell, A. Cooray and 
Rev. D 76, 023507 (200 7) [arXiv:astro-ph/0703375l ; 
B. Jain and P. Zhang, larXiv:0709.2375l [astro-ph]: 
Nesseris and L. Perivolaropoulos, Phys. Rev. D 
103511 (2006) [arXiv:astro-ph/0602053l ; S. Wang, 
Hui, M. May and Z. Haiman, Phys. Rev. 
76, 063503 (2007) arXiv:0705.0l65l [astro-ph]]; 
Tsujikawa, Phys. Rev. D 76, 023514 (2007) 
|arXiv:0705.'T032l [a stro-ph]]; A. F. Heavens, T. D. Kitch- 
ing and L. Verde, arXiv:astro-ph/0703191 ; S. Nesseris 
and L. Perivolaropoulos, Phys . Rev. D 75, 023517 
(2007) [arXiv:astro-ph/0611238] ; J. P. Uzan, Phys . 
Rev. D 59; 123510 (1999) [arXiv:gr-qc/9903004l ; 
V. Sahni and Y. Shtan ov, JCAP 0311, 014 (2003) 

023515 (2005) 



S. 
73 

L. 
D 
S. 



D 71, 



arXiv:astro-ph/0202346 . 
A. Vikman, Phys. Rev. 

[arXiv:astro-ph/0407107l. 

J. Benjamin et al, |arXiv:astro -ph/0703570 C. Shapiro 
and S. Dodelson , Phys. Rev. D 76, 083515 (2007) 
[arXiv :0706.2395 [ast ro-ph]]; L. Amendola, M. Kunz and 
D. Sapone , arXiv:0704.242T] [astro-ph]; A. Refregier et 
!., |arXiv:astro-ph/06100"62| 



[12] 
[13] 

[14] 
[15] 



der, arXiv:0709.1113 [astro-ph]; S. Nesseris and 

L. Perivolaropoulos, larXiv:0710.1092l [astro- ph]; C. Di 

Porto and L. Amendola, arXiv:0707.2686 [astro-ph]; 

Y. Wang, larXiv:0710. 3885l [astro-ph]. 

A. G. Riess et al, |arXiv:astro-ph/0611572| 

T. M. Davis et al, Astrophys. J. 666, 716 (2007) 



[11] A. J. S. Hamilton, |arXiv:astro-ph/ 9708102 E. V. Lin- 



[16 
[17] 
[18] 

[19] 
[20] 
[21] 
[22] 

[23] 



arXiv:astro-ph/0701510 
P. Astier et al., |arXiv:astro-ph/0510447| 
P. Meszaros, Rept. P rog. Phys. 69, 2259 (2006) 
[arXiv:astro-ph /0605208] ; E. W. Liang and B. Zhang, 
Astrophys. J. 633, 611 (2005) [arXiviastro-ph /0504404] ; 
D. Lazzati, G. Ghirlanda, G. Ghisellini, L. Nava, C. Fir- 
mani, B. Morsony and M. C. Begelman, AIP Conf. Proc. 
836 (2006) 513 [arXiv:astro-ph/0602216] ; D. Hooper 
and S. Dodelson, Astr opart. Phys. 27, 113 (2007) 
[arXiv:astro-ph/0512 232 ; F. Y. Wang, Z. G. Dai and 
Z. H. Zhu. larXiv:0706.0938l [astro-ph] . 
M. Bonamente, M. K. Joy, S. 
J. E. Ca rlstrom, E. D. Reese and 
|arXiv:astro-ph/0512349 L. Samushia, G 
B. Ratra. larXiv:0706.1963l [astro-ph] . 
J. P. Uzan, N. Aghanim an d Y. Mellier 
D 70, 083533 (2004) [arXiv:astro-ph/0405620] ; F. De 
Bernardis, E. Giusarma and A. Melchiorri, Int. J. Mod. 
Phys. D 15, 759 (2006) [arXiv:gr-qc/0606029l . 
R. A. Daly, M. P. Mory, C. P. O'Dea, P. Kharb, 
S. Baum, E. J. Guerra and S. G. Djorgovski, 
larXiv:0710.5TT2l [astro-ph]; S. Podariu, R. A. Daly, 
M. P. Mory and B. Ratra, Astrophys. J. 584, 577 (2003) 



J. 
K. 



La Roque, 
S. Dawson, 
Chen and 

Phys. Rev. 



arXiv:astro-ph/0207096 
L. 



018 



72, 



025 



Nesseris and L. Perivolaropoulos, JCAP 0701, 
(2007) [arXiv:astro-ph/0610092] . 
S. Nesseris and L. Perivolaropoulos, Phys. Rev. D 
123519 (2005) [arXiv:astro-ph/0511040] . 
S. Nesseris and L. Perivolaropoulos, JCAP 0702, 
(2007) [arXiv:astro-ph/0612653] . 

W. M. Wood-Vasey et al. [ESSENC E Collaboration], As - 
trophys. J. 666, 6 94 (2007) |arXiv:astro-ph/070104"T| ; 
G. Miknaitis et al. , |arXiv:astro -ph /070 1 043 1 
P. S. Corasaniti and A. Melchiorri. larXiv:0711.4119l 
[astro-ph]; O. Elgaroy and T. Multamaki, Astron. As- 



9 



trophys. 471, 65 (2007) arXiv:astro-ph/0702343 
[24] U. Alam, V Sahni and A. A. Starob insky, JCAP 0702, 

Oil (2007) [arXiv:astro-ph /0612381] . 

[25] Y. Wang and P. Mukherjee, |arXiv:astro-ph/0703780| 
[26] W. J. Percival, S. Cole, D. J. Eisenstein, R. C. Nichol, 

J. A. Peacock, A. C. Pope and A. S. Szalay, 

larXiv:0705.3323l [astro-ph] . 
[27] E. V. Linder, Phys. Rev. D 68, 083503 (2003) 

|arXiv:astro-ph/021230T1 . 
[28] M. Chevalh er and D. Polarski, In t. J. Mod. Phys. D 10, 

213 (2001) [arXiv:gr-qc/0009008l . 
[29] W. Hu, M. Fukugita, M. Zald arriaga and M. Tegmark , 

Astrophys. J. 549, 669 (2001) [arXiv:astro^ph /0006436 ; 

W. J. Percival et al. [The 2dFGRS Team Collabora- 

tion], Mon. Not. Roy. Astron. Soc. 337, 1068 (2002) 



arXiv:astro-ph/0206256 . 
[30] S. Nesseris and L. Perivolaropoulos, Phys. Rev. D 70, 

043531 (2004) [arXiv:astro-ph/0401556T . 
[31] W. H. Press et. al, "Numerical Recipes", Cambridge 

University Press (1994). 
[32] Y. Wang and P. Mukherj ee, Astrophys. J. 650, 1 (2006) 

[arXiv:astro-ph /060405 1 1 . 
[33] D. J. Eisenstein et al. [SDSS Collaboration], Astrophys. 

J. 633, 560 (2005) [arXiv:astro-ph/0501171| . 
[34] B. A. Bass ett and M. Kunz, P hys. Rev. D 69, 101305 

(2004) [arXiv:astro-ph/0312443] . 
[35] A. Mirizzi, G. G. Raffelt and P. D. Serpic o, 

|arXiv:astro-ph/0607415l C. Burrage, larXiv:0711.2966l 



